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We derive and analyze the dynamic equations for polar liquid crystals in two spatial dimensions 
in the framework of classical dynamical density functional theory (DDFT). Translational density 
variations, polarization, and quadrupolar order are used as order-parameter fields. The results 
are critically compared with those obtained using the macroscopic approach of time-dependent 
Ginzburg-Landau (GL) equations for the analogous order-parameter fields. We demonstrate that for 
both the microscopic DDFT and the macroscopic GL approach the resulting dissipative dynamics 
can be derived from a dissipation function. We obtain microscopic expressions for all diagonal 
contributions and for many of the cross-coupling terms emerging from a GL approach. Thus we 
establish a bridge between molecular correlations and macroscopic modeling for the dissipative 
dynamics of polar liquid crystals. 
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I. INTRODUCTION 

Understanding the dynamie processes in liquid erys- 
talline systems from a microseopie point of view is im- 
portant for designing smart materials with novel optieal 
and Theological properties. In fact, the key mechanisms 
of optical displays and switching devices are controlled by 
the dynamic response of liquid crystals to external stim- 
uli [il-Q. Therefore, it is highly relevant to understand 
how these processes change in terms of the molecular in- 
teractions. As first example, the switching dynamics of 
the nematic director upon a change in an external align- 
ment field [J, Q is one of the basic functions in optical 
displays. Second, in microfluidic devices, micromixers [y] 
can be tailored by the response of a liquid crystalline 
system to a rotating external field. 

Nonequilibrium dynamic processes in liquid crystals 
are also interesting from a more fundamental point 
of view since they involve a nontrivial coupling be- 
tween translational and orientational degrees of freedom. 
Therefore, dynamic macroscopic Ginzburg-Landau (GL) 
approaches have been applied to liquid crystalline sys- 
tems in order to obtain the basic dynamic equations on a 
phenomenological level. Apart from direct computer sim- 
ulations of molecular systems [3, Q , much less has been 
done in terms of a microscopic theory. Such a microscopic 
approach which starts from the molecular interactions is 
established by classical density functional theory (DFT) 
[9l-[l3|. DFT can be generalized towards dynamics for 
colloidal particles in solution within the so-called dynam- 
ical density functional theory (DDFT) |14ll6| , which has 
been proven to be a realistic microscopic description for 
many phenomena including liquid crystalline dynamics 

Static DFT constitutes a framework to derive GL equa- 
tions from a microscopic level. The idea is to expand 
the microscopic one-particle density in terms of order- 
parameter fields and to perform a gradient expansion 



[20|-|25[ of the free-energy functional. This leads to a GL- 
like theory, which can also be called a phase-field-crystal 
(PFC) model [2l428J . The phenomenological coupling 
constants of various terms can thereby be expressed in 
terms of generalized moments of molecular correlation 
functions. 

In this paper, we perform a similar analysis for the 
dynamics of liquid crystals by using DDFT, on the one 
hand, and time-dependent GL equations, on the other 
hand. Therefore, we provide a microscopic basis for 
time-dependent GL theory and derive microscopic ex- 
pressions for most of the coupling constants. We do 
this in two spatial dimensions for polar liquid crystals 
by including translational density variations, polariza- 
tion, and quadrupolar order as the basic order-parameter 
fields. We further demonstrate that the resulting dissi- 
pative dynamics can be obtained from a dissipation func- 
tion. This DDFT result holds for constant mobilities as 
well as for mobilities depending on the thermodynamic 
variables. Our work opens the way to study various dy- 
namic processes of two-dimensional polar liquid crystals 
recently observed in experiments |29l - [31| by further nu- 
merical solution of the microscopically justified GL equa- 
tions. 

The paper is organized as follows: in Sec.|IT]we derive 
dynamic equations for polar liquid crystals in two spa- 
tial dimensions from DDFT. A macroscopic derivation 
in terms of time-dependent GL equations and its rela- 
tion to the microscopic dynamics is presented in Sec. IIIII 
We conclude in Sec. IVl 



II. 



MICROSCOPIC DERIVATION OF THE 
DYNAMICS 



Our microscopic derivation of the dynamics uses the 
DDFT equation [13] for collective Brownian motion of 
anisotropic uniaxial colloidal particles and is thus based 



on a static free-energy functional that can be derived 
from static DFT theory. A perturbative functional, 
which uses molecular fluid correlations as input, has re- 
cently been proposed in our previous work [25| for uni- 
axial polar liquid crystalline particles in two spatial di- 
mensions. In the following paragraph, we present this 
functional in a modified form that is more appropriate 
for our further considerations. 



A. Static free-energy functional 

A suitable quantity to describe the state of a system of 
many interacting anisotropic particles is the one-particle 
density p(f,u). It is proportional to the probability den- 
sity to find a particle with orientation ii at position r, 
where u = (ui,M2) = (cos((/3),sin((p)) is the orientational 
unit vector that denotes the orientation of the symme- 
try axis of the respective particle and f= (xi,X2) is the 
center-of-mass position vector. In the following, slightly 
different from Ref. [2^, we choose the parametrization 

p{r,u) = p (l + 4){r) + Pi{r)ui + UtQij{f)uj) (1) 

with the mean number density p and the dimension- 
less order-parameter fields tpi'^^i -Pi(^)i ^^'^ Qiji"^! 
where Einstein's sum convention is used. These order- 
parameter fields are the reduced translational density 



1 



ip{r) = —^ / du(p(f, u) - p) 



(2) 



with iSi denoting the unit circle. The field ipi''^ measures 
translational deviations of p{f, ii) from the mean density 
p. Secondly, the polarization P(r) has the components 



Pi{r) = ^ / dup{r,u)ui 

T^PJsi 



(3) 



and describes the averaged orientation. Finally, the 
traceless and symmetric nematic tensor with the com- 
ponents 



Qij{r) = ^ / dup{r, u) (uiUj - -S.,j) 



(4) 



and the Kronecker delta symbol 5ij describes quadrupo- 
lar ordering. The free-energy functional J-[p{r, u)] is de- 
composed as 

T[p{f, it)] = J"id [p{r, ii)] + -Fcxc [p{r, u)] (5) 

into the ideal rotator gas functional 

/?-Fid[p(f,u)] = fdr fdupir,u){HA^pir,u)) - l) (6) 
Ja Jsi 

with the inverse thermal energy /3 = l/{k-QT), the two- 
dimensional domain A^ and the thermal de Broglie wave- 
length A as well as the the excess free-energy functional 



J^cxc[p(f, u)], which in general is only known approxima- 
tively. Inserting the parametrization ([1} into Eq. ^ and 
performing a Taylor expansion up to fourth order in the 
order-parameter fields yields to the approximation 

/? J-id [p(f, u)] = F,d + Trp /df /id (f) (7) 

Ja 

with the local scaled ideal rotator gas free-energy density 
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(8) 



where 

F,,i = 2npA{HA^p)-l) 
is an irrelevant constant and 

A^ fdr 
Ja 



(9) 



(10) 



is the total area of the domain A. Considering only 
second-order terms in the order-parameter fields and per- 
forming a gradient expansion up to second order, we ob- 
tain from the Ramakrishnan- Yussouff approximation |32l | 



PJ^cKc[pir,u)] = -- dfi dui df2 / dM2 
^ J A JSi J A JSi 

X c'^^\fi -f2,Ml,U2) 

X Ap(fi,ui)Ap(r2,M2) 
with the direct pair-correlation function 

c'^^^(ri,f2,ui,M2) = c'^^\ri ~r2,ui,U2) 



(11) 



(12) 



and the reduced density Ap(r, u) = p{f, u) — p for the 
excess free-energy functional the approximation 



PTcxcipir^u)] 



1 



drfc^dr) 



(13) 



with the local scaled excess free-energy density 

/cxc(r) - Alt// + A2{d,^f + Asidlijf 

+ Biid,^)P, + B2Pr{djQ,j) + B^{d,i^){d,Q,,) 
+ CiPf + C2P^{dlP,) + C^{d,P,f 
+ D^Ql+D2[d,Q,,Y . 

The various coefficients are given by 

^i=8M0(1), 
A2 = -2M0(3), 



(14) 



(15) 
(16) 



^3 = ^M°(5), 

Bi =4(M^i(2) -M?(2)) 
B2 = 2(m}(2)-M^i(2)) 
i33 = -M2_2(3)-M°(3), 
Ci=4Mj(l), 

C3 = -ML2(3) , 

D2 - -M2(3) 
as linear combinations of the moments 



(17) 

(18) 
(19) 
(20) 
(21) 

(22) 

(23) 
(24) 
(25) 



Mr(«)=7rV/di?i?"£f2,(i?) (26) 

Jo 



of the Fourier coefficients 

r2-K r2iT 



,-(2)^m.. -^ 



d^R / d0 c(2) (i?, 0R, 0)e-^('*-+'"^) (27) 







of the direct pair-correlation function c(^)(i?, ^r, 0), for 
which the parametrization 



& '{ri -r2,ui,U2) = c'- 'iR,4>R, 



= .(2) 



(28) 



with f 1 - r 2 = ^u(vr), ""i = ti(Vi) for * = 1,2, (?!)r = 
if — (y9R, and (/) ^ ifi ^ ip2 was used. Equations ([7]) and 
([T3| give a local functional of a polar liquid crystalline 
system reminiscent of a PFC model [19J . 



B. Dynamic equations 

We now derive dynamic equations for ^ = R^ for 
the order-parameter fields ^{r,t), Pi{f^t), and Qij{r,t) 
from dynamical density functional theory. DDFT is con- 
structed to describe the Brownian dynamics of colloidal 
particles in a viscous solvent [3j,[3J] via a time-dependent 
one-particle density field p(f, w,i). This theory was re- 
cently extended to anisotropic Brownian particles with 
orientational degrees of freedom [13, [33, [Sa] ■ It provides 
as a starting point for the case of uniaxial particles in 
two spatial dimensions the DDFT equation [17 1 



dt y dp{r,u,t) ^ 

dip\ ' ' dip Sp{f, u, t) 



(29) 



with the translational short-time diffusion tensor 

Dt(u) = -D||U(8)?i + -D±(l - w®") . (30) 

Here, D|| and D± are the translational diffusion coeffi- 
cients for translation parallel and perpendicular to the 



orientation it, respectively, Dr is the rotational diffusion 
coefficient, ® is the dyadic product, and the symbol 1 
denotes the two-dimensional unit matrix. The two terms 
on the right-hand-side of this DDFT equation for uni- 
axial particles correspond to pure translation and pure 
rotation, respectively. Translational-rotational coupling 
terms, which are especially relevant for screw-like parti- 
cles, do not appear in this DDFT equation, since there is 
no translational-rotational coupling for uniaxial particles. 
Additional terms in the DDFT equation, that regard a 
possible translational-rotational coupling, would have the 
same structure as the present terms, but with only one 
gradient and one angular derivative each instead of two 
gradients or two angular derivatives, respectively [36|. 

Following the analysis of Ref. [13], the functional 
derivative 6J- /5p in the DDFT equation (P5|) has to be ex- 
pressed by functional derivatives of the free-energy func- 
tional with respect to the order-parameter fields ■(/'(r, f). 
Pi (r , t) , and Qij (r , t) , since we parametrized the one- 
particle density p{r,u,t) as well as the free-energy func- 
tional J'[^, Pi, Qij] with these order-parameter fields. In 
the following equations, a large number of functional 
derivatives of the free-energy functional appear. There- 
fore, we shorten the notation by defining the conjugated 
order-parameter fields or thermodynamic forces 



6J- 



^^ = — — with 



s e {p,tlj,Pi,Q,j} 



Using this notation, the equation 
Iirp irp 



np 






(31) 



(32) 



follows by functional differentiation. When performing 
functional derivatives with respect to Qij or Q',, one has 

to notice that Qij as well as Q\, are symmetric and trace- 
less. The interdependence of the elements of these ten- 
sors leads to more complicated derivatives that respect 
the symmetry properties of these tensors. A very useful 
equation in this context is 



5Qki SQ 



kl 



^Q^j 5Q\ 



= 6tk6ji + SjkSti - 5ijS, 



kl 



(33) 



Together with the parametrization ([T} of the one-particle 
density, the relation (^^ can now be inserted into the 
DDFT equation (|29p . The dynamic equations for the 
order-parameter fields are then obtained by an orthogo- 
nal projection that separates the evolution equations for 
the particular order-parameter fields from each other. In 
doing so, the translational density il){r,t) appears to be 
conserved, while Pi(f,t) and Qij(f,t) are not conserved 
due to their association with orientational degrees of free- 
dom. The dynamic equations can thus be written in the 
form 



ij + 8,4 = 

Pr 



•Sf =0 



(34) 
(35) 



Q^, 



<I>g=0 



v- 



(36) field ii G {ip,Pi^Qij} and with the current Jf and the 
quasi-currents ^f and <i>,^. These dissipative currents 



witli ^ = ^ denoting the partial time derivative of the and quasi-currents are given by the expressions 



J7 = 



$ 



(i>'^ = 



- ai(2(l + V)(a.V^) + Qkiid^Qli)) ~ a2Pj{d,P^) 

- a3(2(l + V)(9,Qi-) + ^.(5,/^^) + PMPh + Q^jid,^^)) , 

- 2a,^k{Q^J{^kP^) + Pj{dkQ%)) - a2dk{2{l + 4,){dkPh + Hdui''')) 

- a,(2d,{{l + ^){d,P^)) + 29,((1 + i^MP^)) + d,{P,{d,^^)) + d,{P,{d,i,^)) 

+ 2d,{P,{dkQ%) + QAdkPh)) + «4(2(1 + V^)/^^ + 2P,q5^. - 0,,P,^) , 

- 2aiafe(2(l + ij){dkQ%) + P^idkP)) + P,-(9fe/^^) ~ (5,,Pz(afeP,^) + Q,j{dkiJ^)) 

- Y (49. ((1 + ^){d,iij^)) + H-((i + ^){d,^^)) - 4,5,jai((i + ^)(a,v^)) 

+ Ad,{Pk{d,Pl)) + Ad,{Pk{d,Pl)) ~ ^5,,di{P,{diPl)) 
+ HQH{d,Qli))+d,{QM{d,Qli))-5.,di{Qu,n{diQiJ) 
+ 2^k{Q^,{^lQll)) + '^dk{QM{diQ%)i) 



(37) 



(38) 



PP)+P,P}-5,,PiP}^ 



- J-^ i 



s^,PlPn 



(39) 



where di are the components of the gradient V = {81,82)- 
Four positive coefficients of whom three are independent 
appear in these equations. With the abbreviation A = 
7rp//3, they are defined as 



Ql = 



"3 



D\\+D± 

8A 
DjI - Di_ 

8A 



a2 



04 = 



I?!! + Si:* J. 



8A 



2A 



(40) 



Note that -D|| ^ D± holds for all types of uniaxial parti- 
cles, if the vector ii for the orientation of the symmetry 
axis is chosen properly [37 1. 

At this stage, we emphasize that the DDFT approach 
(|29p a priori contains only three independent mobility 
coefficients, namely the two translational diffusion coef- 
ficients Z?|| and D± and the rotational diffusion coeffi- 
cient Dr. Therefore, all other mobility coefficients for 
the order-parameter fields can be expressed in terms of 
these three basic coefficients. In general, the diffusion co- 
efficients in DDFT are always related to translational or 
orientational degrees of freedom and not to certain order 
parameters, which appear only with the parametrization 
of the one-particle density. The parametrization of the 
one-particle density ([l} in turn does not involve further 
dissipation coefficients. This is in sharp contrast to GL 
theory, where every additional order parameter involves 
at least one new dissipative coefficient - as will be dis- 
cussed in Sec. Illll in more detail. 

Since Eqs. ([57|) - ([5n|) are rather complicated, for nu- 
merical calculations a simpler version of these equations 



I 

might be desirable. Such a simplification is the constant- 
mobility approximation (CMA), where the one-particle 
density in the translational and rotational mobility terms 
of the DDFT equation (|29p is approximated by its mean 
value p: 



|,,.,<,.,,V^ D.,.,V^^Mm 



^pDn 



8^ 6J-[pir,u,t)] 
8(p^ 6p{f,u,t) 



(41) 



With Eq. gT]) instead of the DDFT equation ^, the 
following dissipative currents and quasi-currents are ob- 
tained: 






2ai{8,i;^)^2a^{8jQ%) 



$ 



Q 



^-a2{8M) 



Aa3{8,8jP^) 



^a^P^ , 



2a3 (2(9,5, V-^ ) -^„- (52 V'^)) 



(42) 
(43) 

(44) 



For both the general Eqs. ([37 |) -(|39 | and the much simpler 
constant- mobility Eqs. P^ - (|H| . the explicit forms of the 
conjugated order-parameter fields tp^{f,t), P^{r,t), and 
Q^Jr, t) result directly from the functional derivatives of 
Eqs. ([7]) and ([T^ with respect to the order-parameter 



fields. These functional derivatives are given by 



A b^ 



ir,t) 



QijPj , 



l^<-')-*(« ^ ...... 



P, 



(4 + Ql 



p.p2 



(45) 



(46) 



A oQ 



J-ijri 



U 



^ QijQki 



(47) 



and 



JJ^c: 



~ Bi{d,Pi) - B^{d,djQ,j) , 



Sip 

-, ^-^ CXC , 



-2/3::^(f,i) =Biid,ij) + B2{djQ,j) + 2CiP, 
+ 2C2{dlP,)^2C3{d,djPj), 
B2{d,Pj + djP, - 5,j{diPi)) 



.2/3^1^(r',t) 



B3(2(a.9,V)-%(afV))+4i? 



i^/y 



(48) 
(49) 

(50) 



2D2dk{diQkj + SjQfei - 5ij{diQki)) 



C. Dissipation function 

In the field of linear irreversible thermodynamics [38 
|40| . the dissipative parts of the currents and quasi- 
currents arising in the balance equations for the ther- 
modynamic variables (including, for example, hydrody- 
namic and macroscopic variables) can be derived from 

I 



a dissipation function IH which is quadratic in the ther- 
modynamic forces. Frequently, one uses equivalently the 
entropy production '•JK/T with T denoting the absolute 
temperature |38l . l4l| . The entropy production emerges 
as a source term in the balance equation 



d^jl 






(51) 



for the entropy density a, where j'^ is the entropy cur- 
rent density. The variational derivative of the dissipation 
function with respect to the thermodynamic forces gives 
then the currents and quasi-currents, which are - in lin- 
ear irreversible thermodynamics - by construction linear 
in the thermodynamic forces J38l - l4l| . 

Both, dissipation function and entropy production, are 
maximized close to local thermodynamic equilibrium and 
thus a useful tool in the determination of the dissipa- 
tive currents. This approach has been applied to a large 
number of hydrodynamic and macroscopic systems [38- 
|42| . but is in general not applicable for active systems 
and for systems driven far from equilibrium (compare, 
for example, Refs. [40, 142]) ■ I'^ ^^^i^ more general case, a 
number of additional conditions must be satisfied p3l - |47| 
in order to obtain a Ljapunov functional. Far away from 
equilibrium, the Ljapunov functional is the analogue of 
the dissipation function of linear irreversible thermody- 
namics. 

More precisely, the dissipative currents and quasi- 
currents of the dynamic equations ([M)) - (P5| are given by 



J. 



$f 



$^ 



^ 



5^ 



(S(a,;V'') 



(52) 
(53) 
(54) 



The dissipation function that corresponds to the dissipa- 
tive currents and quasi-currents (|37 l) -(|39 l) of the general 
phase-field-crystal (PFC) model is found to be 



5H(PFC)^ /dr 



ai ((1 + V^)((9,^^)2 + {dkQ%f) + Q^A^k^^){^kQ%) + Q^A^kP^){^kP^) 
+ 2P,{dkP^){dkQl)) + a2 {P.{d,^^){d,P^) + (1 + i^MP^r) 

+ «3((5,:^'^)(iQ.,(a,^^) + P^{d,P^) + P,{d,P^) + 2(1 + mdjQ]^)) 

+ (1 + i^){d,P}f + {d,P}){{l + ^){d,P}) + Q,k{dkP}) + m{dkQ]k)) 

+ \QiAd^Qli){d3Q\i) + \Q^A^^Q'.j){^lQll)) 



(55) 



Together with the dissipative currents and quasi-currents 
(|37|) - ([39)) . this dissipation function constitutes the basic 
rcsuh of this paper. The dissipation function that corre- 
sponds to the currents and quasi-currents (|35])-(|33|) of the 
constant-mobihty approximation (CMA) is much simpler 
and given by 



+ a,{d,P}f + a,{{P}f + 2{Q%f) (56) 



By construction, both dissipation functions ([SS]) and ([55)1 
are positive. This is obvious for Eq. (|56|) . but not mani- 
fest for Eq. ((55)) . 



III. MACROSCOPIC APPROACH: 
GINZBURG-LANDAU DYNAMICS 

In this section, we investigate the Ginzburg-Landau 
(GL) dynamics in the vicinity of the phase transitions 
isotropic to polar ncmatic and isotropic to polar smcc- 
tic. In analogy to the previous section, the GL dy- 
namics is discussed for three types of macroscopic vari- 
ables. These are the smectic density variation p^, which 
is closely related to the complex scalar ip often used to 
describe smectic layering [43, ^^ , the macroscopic polar- 
ization Pi , which becomes important when polar nematic 
and/or polar smectic phases are considered |50l452| . and 
the quadrupolar nematic order parameter Qu, that is 
characteristic of the usual nematic ordering [5j, [SJI • 

We assume that the local formulation of the first law 
of thermodynamics, the Gibbs-Duhem relation, is valid 
I39l44li. It can be written in the form 



TAa = de — /idp — p^dfi,f, — h^ dP,; — QijdSij (57) 

with the absolute temperature T, the entropy density a, 
the energy density e, the chemical potential p, the num- 
ber density p, the chemical potential /z^ associated with 
the layering p^, the thermodynamic force hf associated 
with the macroscopic polarization P,, and the thermo- 
dynamic conjugate Sij of the nematic order parameter 

Throughout the following, we focus entirely on the dis- 
sipative dynamics of the variables associated with the 
additional degrees of ordering, i. e., layering p^, polar 
order Pi, and quadrupolar orientational order Qij. For 
the associated dynamic balance equations, we have one 
dynamic equation each for every hydrodynamic or macro- 
scopic variable. These dynamic equations take the form 
of a conservation law for conserved quantities and are of 
balance equation type for hydrodynamic variables asso- 
ciated with spontaneously broken continuous symmetries 
and for macroscopic variables such as order parameters 



close to a phase transition, 
tions take thus the form 138 



The dynamic balance equa- 
,lil|,li3,[5l[5l[5i 



p^ + dtX, 
P. 



V- 



0, 



yp^o 



Q^ 



z. 



Q 



(58) 
(59) 
(60) 



The currents and quasi-currents X^, K^^, and Z^- are 



introduced via Eqs. ([58|) - (|60|) . Further below, the dis- 
sipative part of their structure will be determined from 
the dissipation function Dl'^'^^K There are no reversible 
currents and quasi-current throughout this paper, since 
flow effects associated with a velocity field v or with a 
density of linear momentum g, are generally not consid- 
ered for the completely overdamped Brownian dynamics 
described by DDFT. We note that the dynamic equation 
associated with the smectic layering is of conservation- 
law type, while the equations for polar and non-polar 
orientational order are balance laws. 

In the spirit of linear irreversible thermodynamics, we 
expand the dissipation function U\^'^^' quadratically in 
the thermodynamic forces p^, hf , and Sij. Those in turn 
have to be determined by taking variational derivatives 



Pi, 
h 



5l_ 

Spij 

^6Pi' 
5F 



5Q^, 



(61) 
(62) 
(63) 



of the suitably chosen generalized potential T with re- 
spect to the variables, where J" has been discussed in 
detail in Refs. [Mill. 

For the dissipation function associated with the three 
types of order considered here, we have to lowest order 
in the gradients 



^0 = <ir(-jijkiSijSki + -aij{diPi,){djPi,) 

+ l^K^i^^ + l^iki{diP4,)Ski 
+ a[j{diPi,)hf + PfhihfSkij ■ 



(64) 



It is usual to consider at first only the lowest-order gra- 
dient terms in the dissipation function. Then, one in- 
spects whether contributions containing more gradients 
are physically relevant. For example, one can always 
add a term containing two more gradients for diagonal 
terms. This was done for the diagonal term ^ 'jijki hi 
Eq. (|M)) , which leads to a relaxation of the nematic or- 
der parameter close to the phase transition [53|, by the 
term ~7yfc/rrm in Eq. (|70[) further below, which contains 
two more gradients and is the dissipative analog of the 
gradient energy of the order parameter Qij . 



From the dissipation function ((64)) . wc obtain for the 
dissipativc currents and quasi-currents the expressions 



XT ^- 



(5fR, 



(GL) 



S{difi^) 



^a,j{djfi^4,) 



~p iP 



Pikl S, 



kl 



Y' 



(5in, 



(GL) 




xP 



5hf 

(GL) 



<5 y- {dj ^i^ ) + hj hj + PikiSki 






V 



{Pkij + Pkji - 5ijl3kii){dk^.^) 

5ij^klmm)Skl 



+ {pL + pL~^^^PmM 



kij 
+ i^khj 



'kji 
- Iklji 



,(65) 
(66) 

(67) 



In a truly isotropic phase, one has only two invari- 
ants: the Kronecker delta 5ij and the totally antisym- 
metric symbol tijk ■ To preserve the symmetries of such 
a system, all the diagonal terms in Eqs. (|M1) - (|571) con- 
tribute, while all off-diagonal coupling terms except for 
one (~a^) vanish. Correspondingly, the property ten- 
sors take the form 



lijki = liSikSji + SjkSii) 



(68) 



aSij , and hij 



a^6,. 



bSij for the 



for the non-vanishing 



(compare also Ref . |57| ) , a 
diagonal terms, and af^ - 
off-diagonal term. 

When comparing this result to Eqs. ([55]) and ([55)) . we 
thus arrive at the conclusion that in a time-dependent 
GL approach we have, even to lowest order in the gradi- 
ents, one diagonal dissipative coefficient for each variable 
entering the dynamics. This has to be contrasted to the 
DDFT approach outlined in the last section, where one 
has only three independent dissipative transport coeffi- 
cients in total (even to higher order in the gradients, com- 
pare the discussion below) . In addition, we find here one 
off-diagonal contribution, ^ a^ , which appears to have 
no analogue in Eqs. ([37| -([39 |) . By direct comparison, we 
find explicitly 7 = 2a4, a = 2ai, and b = 2a4. 

As always, all the dissipative transport coefficients can 
depend on all scalar variables in the system including 
p^ and the temperature T. This general dependence on 
scalar quantities arises partially in the general DDFT re- 
sult ([55]) via the factors {1 + ip) instead of 1, when com- 
pared to the CMA. If one allows also for a dependence 
on vector- and tensor- valued variables, such as the po- 
larization Pi and the quadrupolar order Qij , thus giving 
up the assumption of strict isotropy, the picture outlined 
above changes as follows: the coupling terms between 
the force associated with Sij and the forces associated 
with difi^ and hf can be mediated by the presence of a 
macroscopic polarization Pi. In this case, the property 
tensors fi^ki and /3^ 



3j^f,i take the form 



p^ki = mkPi + suPk) , 

figi^P'^ikkPi+SaPk). 



(69) 



Thus, these dissipative cross-coupling terms can only 
contribute in the presence of a macroscopic polarization. 



Furthermore, they bring along two additional dissipative 
coefficients in a dynamic GL description, while this is not 
the case for the DDFT (compare Eqs. (I32l)-(l39l) of Sec. 
IIIBp . Actually, there appears to be no analogue of the 
contribution ^ Piki in DDFT, while for the contribution 

^ Pikl ^^ f^^^ P^ ^ '^4- 

One can also take into account terms containing more 
gradients in the dissipation function, as it has been done 
in the DDFT approach ([55]) and even for the CMA ((56|) . 
To the next order in the gradients, we obtain 



^(GL) ^ 



•^^ \7^1i]klmn {drnSij ) (9„ Ski ) 

+ ^kMdkhr){dihf) 

+ f3iklmidifJ.^,)idmSkl) 

+ C4(5./^V.)(5fc/if) 

'^ iiklmnK'^nhi ){OmSkl) 
+ Piklmhf {dmSkl) , 



(70) 



In a truly isotropic phase, the following picture emerges 
when Eq. ([70)1 is analyzed: the contribution ~ 'Jijkimn 
is the dissipative analogue of the gradient energy for the 
nematic order parameter. It contains one independent 
material parameter in two spatial dimensions and two 
parameters in three spatial dimensions. In two spatial 
dimensions, we have 7 ~ ai. The tensor bij^i has two 
independent parameters via 



hjki = hSikSji + b2{5ij5ki + Sjk^ii) , 



(71) 



while the tensor Pikim contains one independent param- 
eter 



Piklm = P{SikSlm + SuSkm) ■ 

The same applies to the tensor /3^;,„: 

Kklm = p{SikSlm + SilSkm) ■ 



(72) 



(73) 



All other contributions in Eq. (|70p vanish in a truly 
isotropic phase. Making now an explicit comparison with 
DDFT, we find 61 = 4a3, 62 = a2, and /3 = 03, while 
the contribution ^ /3^ has no analogue in DDFT. The 
same applies to all contributions which simultaneously 
contain only one gradient and are odd in powers of the 
polarization. 

The contributions ~ C,^-/, and 
to contribute as soon as one allows a dependence of the 
property tensors on the polarization Pi. These contri- 
butions can also be associated with the general DDFT 
result given by Eq. ([55]) . We note that all the contribu- 
tions found in Eq. (|55p can also be found in the dynamic 
GL approach when one allows for a dependence of the 
property tensors on the vector- and tensor-valued vari- 



<f,k and -^ifurnn ^^ Eq. GOD Start 



ables. Pi and Qij, used here. However, in contrast to 



DDFT, these dependencies bring along numerous addi- 
tional independent coefBcients. 

Thus, we arrive at the conclusion that the DDFT equa- 
tion ([29)) involves three independent dissipative coefH- 
cients in both the general case and the CMA. The corre- 
sponding terms obtained in the GL framework are associ- 
ated with nine independent coefficients for the analogue 
of the CMA, i. e., for property tensors that do not depend 
on the variables. We also note that there are two cross- 
coupling terms to lowest order in the gradients [see Eq. 
(|64)) ]. which do not exist in the current DDFT picture. 

The overall picture that emerges is therefore the fol- 
lowing. In a dynamic GL approach, there is at least 
one independent dissipative coefficient for every dissipa- 
tion channel (every order-parameter field) entering the 
description. In addition, one finds frequently dissipa- 
tive cross-coupling terms that bring along further co- 
efficients - in particular, if one considers a dependence 
of the dissipative property tensors on the macroscopic 
variables. This can be contrasted to the current DDFT 
picture, where one has only two dissipation channels of 
diagonal nature, namely translational and rotational dif- 
fusion. In the present version of DDFT, there are also no 
independent dissipative cross-coupling terms. These ob- 
servations clearly call for a generalization of the current 
DDFT equation to incorporate processes that allow for 
additional dissipation on a microscopic level. 



IV. CONCLUSIONS AND POSSIBLE 
EXTENSIONS 

In conclusion, we have proposed both microscopic and 
macroscopic theoretical descriptions for the dynamics of 
polar liquid crystals in two spatial dimensions. The mi- 
croscopic theory is derived from DDFT, while the macro- 
scopic formulation is based on time-dependent GL the- 
ory. We have done this by including translational den- 
sity variations, polarization, and quadrupolar order as 
the basic order-parameter fields. Most but not all phe- 
nomenologically possible couplings of GL theory occur 
also in the DDFT approach. These couplings are derived 
from a microscopic approach and the associated coupling 
parameters can be expressed as generalized moments of a 



molecular correlation function. We further demonstrated 
that the whole dynamics can be obtained from a dissipa- 
tion function. 

Our theoretical framework can be used for a further 
exploration of various dynamic processes of polar liquid 
crystals. This requires numerical solutions of the mi- 
croscopically justified GL equations following numerical 
schemes proposed earlier [581 i5j| . 

For future work, it is challenging to construct a gen- 
eralized DDFT which explicitly contains the momentum 
field as appropriate for molecular dynamics or systems 
in flow fields J6C| . This turns out to be much more diffi- 
cult than the traditional DDFT approach for simple over- 
damped Brownian dynamics. But in principle the way of 
generalization was explored by Tarazona, Marconi and 
Melchionna |6ll - [63| and by Archer [6J, i65| for molecular 
dynamics. An alternative derivation is based on projec- 
tor techniques [66| leading to a hydrodynamical density 
functional theory |16l . |67[. Additional dynamic expres- 
sions for a colloidal liquid under shear flow were recently 
discussed in Ref. [60] ■ Furthermore, a phase- field-crystal 
model coupled to fiow was considered by Voigt and co- 
workers 631 (see also Ref. [631 )■ In these extensions, 
one will presumably obtain non-vanishing microscopic 
expressions for phenomenological dynamic terms caused 
by the existence of a momentum or velocity field. These 
come macroscopically mainly in two groups, namely con- 
tributions leading to a flow alignment associated with re- 
versible currents coupling extensional flow (symmetrized 
velocity gradients) to orientational degrees of freedom 
pll 15% [Sa . [70l - [76| and coupling terms between exten- 
sional flow and variations of the moduli of nematic, smec- 
tic, and columnar order [H [H [H El, [z3, Izl • 
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